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APPARATUS FOR AND METHOD OF MEASURING JITTER 



BACKGROUND OF THE INVENTION 

The present invention relates to an apparatus for and a method of 
5 measuring period jitter that are applied to a measurement of jitter of, for 
example, a microprocessor clock. 

A time interval analyzer and/or an oscilloscope have conventionally 
been used for the measurement of period jitter. The method of these 
apparatus is called Zero-crossing Method, in which, as shown in Fig. 1, a 

10 clock signal (a signal under measurement) x(t) from, for example, a PLL 
(Phase-Locked Loop) under test 11 is supplied to a time interval analyzer 
12. Regarding a signal under measurement x(t), a next rising edge 
following one rising edge fluctuates against the preceding rising edge as 
indicated by dotted lines. That is, a time interval T p between the two rising 

15 edges, namely a period fluctuates. In the Zero-crossing Method, a time 
interval between zero-crossings (period) of the signal under measurement is 
measured, a fluctuation of period is measured by a histogram analysis, and 
its histogram is displayed as shown in Fig. 2. A time interval analyzer is 
described in, for example, "Phase Digitizing Sharpens Timing 

20 Measurements" by D.Chu, IEEE Spectrum, pp. 28-32, 1988, and "A Method 
of Serial Data Jitter Analysis Using One- Shot Time Interval Measurements" 
by J.Wilstrup, Proceedings of IEEE International Test Conference, pp. 
819-823, 1998. 

In addition, Tektronix, Inc. and LeCroy Co. have recently been 
25 providing digital oscilloscopes each being able to measure a jitter using an 
interpolation method. In this jitter measurement method using the 
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interpolation method (interpolation-based jitter measurement method), an 
interval between data having signal values close to a zero-crossing out of 
measured data of a sampled signal under measurement is interpolated to 
estimate a timing of zero-crossing. That is, in order to measure a 
5 fluctuation of period, a time interval between zero-crossings (period) is 
estimated using a data interpolation with a small error. 

That is, as shown in Fig. 3, a signal under measurement x(t) from the 
PLL under test 11 is inputted to a digital oscilloscope 14. In the digital 
oscilloscope 14, as shown in Fig. 4, the inputted signal under measurement 

10 x(t) is converted into a digital data sequence by an analog-to-digital 

converter 15. A data- interpolation is applied to an interval between data 
having signal values close to a zero-crossing in the digital data sequence by 
an interpolator 16. With respect to the data-interpolated digital data 
sequence, a time interval between zero-crossings is measured by a period 

15 estimator 17. A histogram of the measured values is displayed on a 
histogram estimating part 18, and a root-mean-square value and a 
peak-to-peak value of fluctuations of the measured time intervals are 
obtained by an RMS & Peak-to-Peak Detector 19. For example, in the case 
in which a signal under measurement x(t) has a waveform shown in Fig. 5A, 

20 its period jitters are measured as shown in Fig. 5B. 

In the jitter measurement method by the time interval analyzer 
method, a time interval between zero-crossings is measured. Therefore a 
correct measurement can be performed. However, since th e re is, in this 
jitter measurement m e thod, a dead time wh e n no m e asurem e nt can bo 

25 performed after one period mea s urement i s p e rformed. H owever, because 
this method repeatedly measures jitter but includes an intermediate 




-3- 

dead-time between measurements, t here is a problem that it takes a long 
time to acquire a number of data that are required for a histogram analysis. 
In addition, in an interpolation-based jitter measurement method in which a 
wide-band oscilloscope and an interpolation method are combined, there is a 
5 problem that a histogram of jitter cannot accurately be estimated, and a jitter 
values are overestimated (overestimation). For example, for a 400 MHz 
clock signal the time interval analyzer method measure a root-mean-square 
value of jitter as 7.72 ps while the interpolation method measures, a 
root-mean-square of 8.47 ps, that is larger than the value estimated by the 

10 time interval analyzer method. 

On the other hand, inventors of the present invention have proposed 
a method of measuring a jitter as described below in an article entitled 
"Extraction of Peak-to-Peak and RMS Sinusoidal Jitter Using an Analytic 
Signal Method" by T.J.Yamaguchi, M.Soma, M.Ishida, and T.Ohmi, 

15 Proceedings of 18th IEEE VLSI Test Symposium, pp.395-402, 2000. That 
is, as shown in Fig. 6, an analog clock waveform from a PLL (Phase locked 
loop) circuit under test 11 is converted into a digital clock signal x c (t) by an 
analog-to-digital converter 22, and the digital clock signal x c (t) is supplied 
to a Hilbert pair generator 24 acting as an analytic signal transforming part 

20 23, where the digital clock signal x c (t) is transformed into an analytic signal 

Zc(t). 

Now, a clock signal Xc(t) is defined as follows. 
x c (t)=Ac cos(27tf c t+e c -A(|>(t)) 

25 

The Ac and the f c are nominal values of amplitude and frequency of a clock 



signal respectively, the 9 C is an initial phase angle, and the A<|>(t) is a phase 
fluctuation that is called an instantaneous phase noise. 

Signal components around a fundamental frequency of the clock 
signal x c (t) are extracted by a bandpass filter (not shown) and are 
Hilbert-transformed by a Hilbert transformer 25 in the Hilbert pair generator 
24 to obtain the following equation. 

X c (t) = H[X C (t)] = A c sin(2icf c t + 9 C - A<|>(t)) 

Then, an analytic signal z c (t) having x c (t) and x c (t) as a real part and an 
imaginary part, respectively, is obtained as follows. 

-Z 6 (t)=* e(t)+X c (/) 

~ A c cos(2^ c f I 6 C &<j>(t) I jA c mn(2nf c t I B c A flr)) 

Zr (i)=Xr(t)+\ X c (t) 

= A c cos(27tf c t + 0 C - Aq>(t)) + jA c sin(27if c t + 9 C - A(p(t)) 

From this analytic signal z c (t), an instantaneous phase 0(t) of the 
clock signal Xc(t) can be estimated by the instantaneous phase estimator 26 
as follows. 

0(t)=[27if c t+e c -A<|)(t)] mod 2tt [rad] 

A linear phase is removed from this instantaneous phase 0(t) by a linear 
phase remover 27 to obtain a phase noise waveform A<|)(t). That is, in the 
linear phase remover 27, a continuous phase converting part 28 applies a 
phase unwrap method to the instantaneous phase 0(t) to obtain a continuous 
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instantaneous phase 0(t) as follows. 

0(t)= 27tf c t+0 c -A(|)(t) [rad] 

5 The phase unwrap method is shown in "A New Phase Unwrapping 
Algorithm" by Jose M. Tribolet, IEEE Trans. Acoust., Speech, Signal 
Processing, vol. ASSP-25, pp. 170-177, 1977 and in "On Frequency-Domain 
and Time-Domain Phase Unwrapping" by Kuno P.Zimmermann, Proc. IEEE, 
vol. 75, pp. 519-520, 1987. 

10 An instantaneous linear phase of a continuous instantaneous phase 

0(t), i.e., a linear instantaneous phase [27tf c t+0 c ] of a jitter-free ideal signal 
is estimated by a linear phase evaluator 29 using a linear trend estimating 
method. That is, an instantaneous linear phase of a continuous 
instantaneous phase 0(t) is estimated by applying a linear line fitting by least 

15 squares method to the above continuous phase 0(t). This estimated linear 
phase [27if c t+0 c ] is subtracted from the continuous phase 0(t) by a 
subtracting part 31 to obtain a variable term A<|>(t) of the instantaneous phase 
@(t) , i.e., an instantaneous phase noise waveform as follows. 

20 9(t)=A<|>(t) 

The instantaneous phase noise waveform A<|>(t) thus obtained is inputted, 
after having been sampled by the zero-crossing sampler 34, to a 
peak-to-peak detector 32 as a timing jitter sequence A(f>[n], where a 
25 difference between the maximum peak value max (A(|>[k]) and the minimum 
peak value min (A(|)[k]) of the A(|)[n] (=A(|)(nT)) is calculated to obtain a peak 



value (peak-to-peak value) A§ pp of timing jitter as follows. 



A(j) pp : 



,= max (A<|>[k])- min (A(j)[k]) 

k k 



In addition, the in s tantaneous t iming jitter sequence A<J)[n] is also inputted to 
a root-mean-square detector 33, where a root-mean-square (RMS) value of 
the timing jitter sequence A<))[n] is calculated using following equation to 
obtain a root-mean-square value A^rms of timing jitters. 



This method is referred to as the A<|> method, since a peak value of 
timing jitter (peak-to-peak value) and a root-mean-square value of timing 
jitters are obtained from the instantaneous phase noise waveform A<J)(t) . 
Further, an instantaneous phase noise waveform A<J)(t) is sometimes written 
as a instantaneous phase noise A(|>(t) or a phase noise waveform A(|>(t). 

According to the A<|> method, a timing jitter can be measured at high 
speed with relatively high accuracy. 

It is an object of the present invention to provide an apparatus for 
and a method of measuring a jitter that can measure a period jitter in a short 
period of time and with high accuracy, namely an apparatus for and a 
method of measuring a jitter that can measure jitter values compatible with 
those measured by the conventional time interval analyzer method. 
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SUMMARY OF THE INVENTION 

According to the present invention, an instantaneous phase noise 
waveform of a signal under measurement can be obtained, the instantaneous 
phase noise waveform is sampled at a timing close to each zero-crossing 
5 point (approximated zero-crossing point) of the signal under measurement to 
estimate a timing jitter sequence of the signal under measurement, a 
difference sequence of this timing jitter sequence is calculated to measure a 
period jitter sequence, and values of this period jitter sequence is corrected 
by multipling a ratio of a fundamental period of the signal under 
10 measurement and time interval values between the approximated 
zero-crossing points. 

A principle of the present invention will be explained below. An 
analytic signal of a fundamental cosine wave x(t) of an input signal (may 
sometimes be written as a signal under measurement) is given by an 
15 equation (1) as follows. 

z(t) = x(t)+jH[x(t)] 

= A cos(2nf 0 t + 6 - A^(0) + jA sin(2nf 0 t + 6 - A#/)) ^ ^ 



In this case, f 0 is a fundamental frequency of the signal under measurement, 
and is equal to fo=l/T 0 7 .where T o is a fundamental period of the signal 
20 under measurement. An instantaneous frequency of z(t) is expressed by 
the following equation. 

1 = <p(0 = 1 _x(t) H'[x(t)]-x\t)H[x(t)] 
T 0 +J In In 
1 ( T \ 
T 0 \ 2n J 



x 2 (t) + H 2 [x(t)] 



[Hz] (2) 
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Therefore, the following equation is given. 




[sec] (3) 



5 



A timing jitter sequence is obtained by sampling an instantaneous phase 
noise waveform A<J>(t) at a timing (referred to as an approximated 
zero-crossing point) closest to each zero-crossing point of a real part x(t) of 
an analytic signal z(t), and then it is assumed that the sampling interval 
10 Tk,k+i of the approximated zero-crossing points is equal to a fundamental 
period T 0 . In this case, a period jitter J is obtained, as shown by the 
following equation, as a difference sequence of a timing jitter sequence. 



In the equation (4), the division of 27t/To is performed to convert the 
unit of radian into unit of second. 

However, s ince the in s tantaneou s phas e nois e waveform A4>(t) is 
sampled by approximating a zero crossing point of th e signal under 
20 measurement to a sampling point closest to each zero crossing point, as 
shown in Fig. 7, an approximat e d z e ro cro s sing point x at which tho 
instantan e ous phas e noi se wav e form A$(t) is sampled is s hifted from a 
corresponding z e ro cro s sing point of th e signal und e r measur e m e nt o. 



In 



[sec] 



(4) 



15 
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That i s , an interval betwe e n approximat e d zero crossing points x is different 
from th e fundam e ntal p e riod T or Since the instantaneous phase noise 
waveform is sampled at sampling points closest to ideal zero-crossing 
points, as shown in Fig. 7. an approximated zero-crossing point x at which 
5 the instantaneous phase noise waveform A$(t) is sampled is deviated from a 
corresponding ideal zero-crossing point of the signal under measurement o. 
That is. a time interval between approximated zero-crossing points x is 
different from the fundamental period Tn. Therefore, when the equation M) 
is calculated based on th e timing jitter sequ e nc e obtained by sampling the 

10 instantaneous phas e nois e wav e form A(|)(t) at thos e approximated 

zero - crossing points, a period jitter cannot be obtained with high accuracy. — 
Particularly in the case where a sampling period i s large and an 
ov e r sampling rate i s small, a measur e m e nt e rror of p e riod jitt e r b e come s 
larg e . — In order to d e cr e as e this e rror, data of mor e than 10 point s 

15 (ov e r sampling rate is 5) or so p e r period Tqare r e quir e d. Therefore, when 
the equation (4) is used, a period jitter cannot be estimated with high 
accuracy. Particularly in the case where a sampling period is large and an 
over sampling ratio is small, a measurement error of period iitter becomes 
large. In order to achieve a smaller error of this type, data of more than 10 

20 points per period Tn (over-sampling ratio is 5) or so are required. 

According to the present invention, a correction is perform e d, as 
shown in an e quation (5), by multiplying the e quation (1) by the ratio of the 
fundamental period To to the int e rval a correction can be realized using 
equation (5^ bv multiplying equation (4) bv the ratio of the fundamental 

25 period T n to the time interval b etween the approximated zero-crossing points 

Tk,k+1- 
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In 



I T 



[sec] 



(5) 



In this cas e , To A^fie-^ ^ repres e nts a corr e ctive term for th e diff e r e nc e bas e d 
5 approximation ( e quation (4)) of th e instantan e ous phas e noise differential. — 
A period jitter can be obtained with high accuracy u s ing this corrective term. 
As s hown in Fig. 7, the time interval b e tween approximated 
zero - crossing points can be obtained by obtaining a timing sequence t[k] of 
the approximated zero - crossing points and then by calculating their 
difference s . where the term. To/ Tv v+i corrects the instantaneous phase noise 
difference using the difference-based approximation (equation (4)). A 
period jitter can be obtained with high accuracy using this correction term. 
As shown in Fig. 7. the time interval TV kjH between approximated 
zero-crossing points can be obtained by differentiating a timing sequence 
15 t[k] at the approximated zero-crossing points. 



10 



T k)k+ i=t[k+l]-t[k] 



[sec] 



(6) 



In addition, the fundamental period T 0 of the signal under 
20 measurement may be obtained from an inclination 2tu/T 0 of the linear 
instantaneous phase or may be directly obtained from the signal under 
measurement. 

When the J[k] is corrected by the correction term as shown in Fig. 8, 
an estimation error (shown in parentheses) between an estimated jitter value 
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and an ideal value with respect to either of an RMS value Jrms and a 
peak-to-peak value Jpp of period jitter can be reduced. The data shown in 
Fig. 8 were obtained through a computer simulation applied to a signal 
having sinusoidale jitter. Particularly, when the over-sampling ratio is 
5 small, its effect becomes large. In Fig. 8, "A(|> method" shows a case where 
a correction by To/Tk, k +i is not applied, and " Corrected A§ method" shows 
a case where the correction by T 0 / Tk, k+i is applied. Fig. 9 shows an 
experimental result when a real waveform is used. Fig. 9 shows measured 
jitter values for various numbers of sampling points per period T 0 . In the 

10 case of the A<|) method (indicated by x) where the correction by T 0 /Tk, k +i is 
not applied, a peak-to-peak value of period jitter is especially over estimated 
as the over-sampling ratio is getting smaller. On the other hand, according 
to the present invention, a peak-to-peak value of period jitter can be 
obtained correctly, by using the above correction term. Particularly, when 

15 the over-sampling ratio is small, its effect becomes large. For example, in 
the example shown in Fig. 9B, the error can be corrected by 8% in the case 
of 8 points per period (over-sampling ratio is 4), and by 18% in the case of 3 
points per period (over-sampling ratio is 1.5). As a result, according to the 
present invention, it has become possible to obtain a period jitter up to the 

20 case where the over-sampling ratio is 1 .5 using the A(|) method. This means 
that, if the sampling period is the same, a jitter of a signal under 
measurement having higher frequency can be measured more accurately. 

In addition, a period for obtaining a period jitter may be m periods 
(m=0.5, 1, 2, 3, • •). That is, by taking the m=0.5 period, a difference 

25 between two timing jitter values at a rising (or falling) zero-crossing point 
and a next falling (or rising) zero-crossing point may be obtained. Also, by 
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taking the m=2 periods, a difference between two timing jitter values at a 
rising (or falling) zero-crossing point and a rising (or falling) zero-crossing 
point after two rising (or falling) points from the first rising (or falling) 
zero-crossing point may be obtained. By calculating a root-mean-square 
5 and a difference between the maximum value and the minimum value of the 
period jitter data measured in this manner, an RMS value Jrms and a 
peak-to-peak value Jpp of the period jitter can be obtained by the following 
equations , respectively. 

10 

Mm 

Jrms = J— U 2 M [sec] (7) 

VMk=i 



J = max( J[k]) - mm(J[k]) [sec] (8) 

rr k k 

15 

In this case, M is the number of sampl e s of th e measur e d period jitt e r data. 
Fig. 10 show s a histogram (Fig. 10B) of the p e riod jitter measured by the 
m e thod and a histogram (Fig. 10A) m e a s ur e d by tho conventional time 
interval analyzer so that the hi s togram of th e A§ m e thod can be compared 
20 with the histogram of the conventional tim e interval analyzer. — In addition, 
Fig. 11 shows an RMS value and a p e ak to p e ak value of th e period jitt e r 
m e asured by th e Corr e cted A(|) m e thod so that th o s e valu e s of th e Corr e cted 
A(|) m e thod can b e compared with th e corresponding values m e asur e d by tho 
tim e int e rval analyzer, respectiv e ly. — H e r e , th e peak to peak valu e Jpp of tho 
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ob se rv e d p e riod jitter is substantially proportional to a square root of 
logarithm of th e numb e r of e v e nts (the number of zero crossings). In this 
case. M is the number of samples of the measured period jitter sequence. 
Fig. 10 shows a histogram (Fig. 10B) of the period jitter measured by the 
5 corrected A6 method and a histogram (Tig. 10 A) measured by the 

conventional time interval analyzer. Thus the histogram estimated bv the 
corrected Ad) method can be compared with the histogram of the 
conventional time interval analyzer. In addition. Fig. 11 shows an RMS 
value and a peak-to-peak value of the period jitter measured by the corrected 

10 A([) method. The corrected A(|) method provides both measured 

peak-to-peak and RMS jitter values that are comparable to the time interval 
analyzer method. Note that the peak-to-peak value Jp p of the observed 
period jitter is substantially proportional to a square root of logarithm of the 
number of events (the number of zero-crossings V F or example, in the case 

15 of approximately 5000 events, Jpp =45 ps is a correct value. A J PP error in 
Fig. 11 is shown assuming that 45 ps is the correct value. As shown in 
Figs. 10A, 10B, and 11, the Corrected A(|) method can provide measured 
jitter values compatible with the conventional measurement method. 

Furthermore, the A§ method can measure a cycle-to-cycle period 

20 jitter and a period jitter simultaneously. A cycle-to-cycle period jitter Jcc 
is a period fluctuation between adjacent clock cycles, and is expressed by an 
equation (9). 



25 



J cc [k]=T[k+l]-T[k] 

=(T 0 +J[k+1])-(T 0 +J[k]) 
=J[k+l]-J[k] [sec] 
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Therefore, by obtaining differences between the period jitter data measured 
as described above and then calculating their root-mean-square and a 
difference between the maximum value and the minimum value, an RMS 
5 value Jcc,rms and a peak-to-peak value Jccjpp of cycle-to-cycle period jitter 
can be obtained by equations (10) and (11), respectively. 

Jcc,RMS = l\tjccW [sec] (10) 

10 ^cc ) pp= m ^(^ccTO- m p(' / cc[^]) [ sec ] (11) 

k k 

In this case, L is the number of samples of the measured cycle-to-cycle 
period jitter data. 

1 5 BRIEF DESCRIPTION OF THE DRAWINGS 

Fig. 1 is a diagram showing a period jitter measurement using a 
conventional time interval analyzer; 

Fig. 2 is a diagram showing a histogram of the measured values; 
Fig. 3 is a diagram showing a jitter measurement using a 
20 conventional digital oscilloscope; 

Fig. 4 is a diagram showing a configuration of a jitter measurement 
part in Fig. 3; 

Fig. 5A is a diagram showing a waveform of a signal under 
measurement; 
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Fig. 5B is a diagram showing measured period jitter of the signal 
under measurement; 

Fig. 6 is a diagram showing a functional configuration of a jitter 
measurement apparatus based on the A<|> method previously proposed by the 
5 inventors of the present invention; 

Fig. 7 is a diagram showing a discrepancy between a zero-crossing 
point and an approximated zero-crossing point; 

Fig. 8 is a diagram showing a comparison of measured 
root-mean-square values and peak-to-peak values of period jitter between 
10 the A<|> method and the method by which the period jitter is corrected (the 
method of the present invention); 

Fig. 9A is a diagram showing an example of the effect of a correction 
term for an RMS period jitter estimation; 

Fig. 9B is a diagram showing an example of the effect of a correction 
1 5 term for a peak-to-peak period jitter estimation; 

Fig. 10A is a diagram showing a histogram of period jitter measured 
by the conventional time interval analyzer method; 

Fig. 10B is a diagram showing a histogram of period jitter measured 
by the Corrected A<|> method according to the present invention; 
20 Fig. 11 is a diagram showing RMS values and peak-to-peak values of 

period jitter measured by the time interval analyzer method and the 
Corrected A§ method; 

Fig. 12 is a block diagram showing a functional configuration of an 
embodiment of the present invention; 
25 Fig. 13 is a diagram showing a functional configuration of a partially 

modified example of the present invention; 
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Fig. 14 is a diagram showing another specific functional 
configuration of the analytic signal transforming part 23; 

Fig. 15 is a flow-chart showing a procedure of an embodiment of a 
method according to the present invention; 
5 Fig. 16 is a block diagram showing a portion of another embodiment 

of the apparatus according to the present invention; and 

Fig. 17 is a flow-chart showing a portion of another embodiment of a 
method according to the present invention. 

10 DETAILED DESCRIPTION OF THE EMBODIMENT 

Fig. 12 shows an embodiment of the present invention. Portions in 
Fig. 12 corresponding to those in Fig. 6 have the same reference numbers 
affixed thereto as those in Fig. 6, and duplicated explanations for those 
portions will be omitted. According to the present invention, a signal 

15 under measurement x c (t) is inputted to a timing jitter estimator 39, and a 
timing jitter sequence of the inputted signal under measurement is obtained. 
In this embodiment, an instantaneous phase noise waveform A(|>(t) obtained 
from a linear phase remover 27 is sampled by a zero-crossing sampler 43 at 
timings closest to zero-crossing points of a real part Xc(t) of an analytic 

20 signal z c (t) to obtain a timing jitter sequence. For this reason, a real part 
x c (t) of an analytic signal from an analytic signal transforming part 23 is 
inputted to a zero-crossing point detecting part 45. 

Next, approximated zero-crossing points are detected by the 
zero-crossing point detecting part 45. That is, the maximum value of a 

25 waveform of an inputted real part x c (t) is defined as 100% level, and the 

minimum value is defined as 0% level to calculate 50% level V (50%) of the 
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difference between the 100% level and the 0% level as a zero-crossing level. 
A difference between a sample value and 50% level V (50%) and a 
difference between its adjacent sample value and 50 % level V (50%), i.e., 
(x c (j-l)-V (50%)) and ( x c (j)-V (50%)) are calculated, and further a product 
5 of those difference values (x c (j-l)-V(50%)) x ( x c (j)-V (50%)) is calculated. 
When the Xc(t) crosses 50% level, i.e., zero level, the sign of its sample 
value x c (j-l) -V(50%) or Xc(j) -V(50%) changes from a negative sign to a 
positive sign or from a positive sign to a negative sign. Therefore, when 
the product is negative, it is detected that the x c (t) has passed the zero level, 

10 and a time point j-1 or j at which a smaller absolute value of the sample 

value x B (i 1) or x B (i) is d e tect e d is obtained as an approximate v alue Xcfi-l) - 
V(50%) or Xcti) - Vf 50%) is detected: an approximate zero crossing point. 
At each approximated zero-crossing point, a sampling pulse is supplied to 
the zero-crossing sampler 43. 

15 A sample value sequence outputted from the zero-crossing sampler 

43, namely a timing jitter sequence, is inputted to a difference calculating 
part 46, where a difference sequence of the timing jitter sequence is 
calculated. That is, the equation (4) is calculated with respect to the 
inputted A<|>[k] and A<|>[k+1], and upon every update of k, the equation (4) is 

20 calculated to estimate a period jitter sequence. In addition, a timing 
sequence t[k] at each sampling time in the zero-crossing sampler 43 is 
inputted to a zero-crossing interval calculating part 47, where a time interval 
Tk,k+i between approximated zero-crossing points is obtained through a 
calculation of the equation (6). 

25 In addition, an instantaneous linear phase from the linear phase 

remover 27, namely a linear phase component from the linear phase 
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estimating part 29 in Fig. 6 is inputted to a fundamental period estimating 
part 48, where a fundamental period To is obtained from the inclination 
2tt/To of the instantaneous linear phase. This fundamental period To may 
be obtained from the signal under measurement itself by inputting the signal 
5 under measurement from the AD converter 22 to the fundamental period 
estimating part 48. Alternatively, when the fundamental period To of the 
signal under measurement is known in advance, that value T 0 may be stored 
in advance in the fundamental period estimating part 48. 

A period jitter sequence from a difference calculating part 46, an 

10 approximated zero-crossing point interval T^k+i from the zero-crossing 
interval calculating part 47, and a fundamental period T 0 from the 
fundamental period estimating part 48 are inputted to a corrector part 49, 
where each period jitter in the period jitter sequence is multiplied by 
To/Tk,k+i, namely the equation (5) is calculated to obtain a corrected period 

15 jitter sequence. 

This corrected period jitter sequence is directly supplied to a 
cycle-to-cycle period jitter estimating part 52, and at the same time, a period 
jitter sequence that is the corrected period jitter sequence delayed by its one 
element (one period jitter) is also supplied to the cycle-to-cycle period jitter 

20 estimating part 52. 

The cycle-to-cycle period jitter estimating part (differentiator) 52 
calculates a difference sequence of the period jitter sequence at each time 
point k using the equation (9) to obtain a cycle-to-cycle period jitter 
sequence. 

25 This embodiment is a case where the corrected period jitter sequence 

from the correcting part 49 and the cycle-to-cycle jitter period jitter 
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sequence are switched by a switch 53 so that one of those corrected period 
jitter sequence and the cycle-to-cycle jitter period jitter sequence can 
selectively be supplied to a jitter detecting part 54. 

In this case, there are provided in the jitter detecting part 54 a 
5 peak-to-peak detecting part 32 for obtaining a difference between the 
maximum value and the minimum value of the inputted jitter sequence, an 
RMS detecting part 33 for calculating a root-mean-square (RMS) value of 
the inputted jitter sequence, and a histogram estimating part 18 for obtaining 
a histogram of the inputted jitter sequence. 

10 In the state that the switch 53 is connected to the output side of the 

corrector part 49, the corrected period jitter sequence is inputted to the jitter 
detecting part 54. In this case, the equation (8) is calculated with respect to 
the period jitter sequence by the peak-to-peak detecting part 32 to obtain a 
peak-to-peak value J pp of period jitter, the equation (7) is calculated by the 

15 RMS detecting part 33 to obtain an RMS value Jrms of period jitter, and a 
histogram of period jitter is obtained by the histogram estimating part 18, 
and then those obtained values are outputted to be displayed on, for example, 
a display part (not shown). 

In the state that the switch 53 is connected to the output side of the 

20 p e ak to p e ak p e riod jitt e r e stimating part 52, the p e ak to peak period jittor 
cvcle-to-cvcle period jitter estimating part 52. the cvcle-to-cvcle period jitter 
sequence is inputted to the jitter detecting part 54. In this case, the 
equation (11) is calculated by the peak-to-peak detecting part 32 to obtain a 
peak-to-peak value J C c,p P of cycle-to-cycle period jitter, and the equation (10) 

25 is calculated by the RMS detecting part 33 to obtain an RMS value Jcc.rms 
of cycle-to-cycle period jitter, and further, a histogram of cycle-to-cycle 
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period jitter is estimated by the histogram estimating part 18. Then those 
obtained values are also outputted and are displayed on the display part if 
necessary. 

In the configuration shown in Fig. 12, the cycle-to-cycle period jitter 
5 estimating part 52 and the switch 53 may be omitted to supply the corrected 
period jitter sequence from the correcting part 49 directly to the jitter 
detecting part 54. In addition, the switch 53 may be omitted to supply the 
cycle-to-cycle period jitter sequence from the cycle-to-cycle period jitter 
estimating part 52 directly to the jitter detecting part 54. Furthermore, the 

10 jitter detecting part 54 may only include any one or two of the peak-to-peak 
detecting part 32, the RMS detecting part 33, and the histogram estimating 
part 18. In order to detect approximated zero-crossing points for sampling 
the instantaneous phase noise A<|>(t), not only the real part signal is utilized 
but also the signal under measurement itself or its fundamental wave 

15 component may be utilized. 

As indicated by dashed lines in Fig. 12, a clock signal from the PLL 
under test 1 1 can be supplied to the AD converter 22 via a waveform clipper 
56 to make amplitude of the clock signal constant. In this arrangement, a 
jitter can corr e ctly accurately b e measured in the state that the phase noise 

20 waveform A<|>(t) is not influenced by an-amplitude modulation component . 
components . This process for making amplitude of the input signal 
constant removing amplitude modulation components from the input signal 
may be performed at the output side of the AD converter 22. 

As shown in Fig. 13, a process of an analytic signal transforming 

25 part 23 for transforming an input signal to an analytic signal z c (t) comprises 
the steps of transforming a digital input signal from the AD converter 22 
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into a both-sided spectrum signal in frequency domain using, for example, 
Fast Fourier Transform (FFT) by a tim e domain to fr e qu e ncy domain 
transforming part 61, cutting off n e gativ e fr e qu e ncy compon e nts from th e 
both sided spectrum signal by a bandpass filter 62 to extract, at the s ame 
5 time, frequency components around a po s itive fundam e ntal fr e qu e ncy of the 
input clock signal by the time domain to the frequency domain transforming 
part 6L making negative frequency components of the both-sided spectra 
zeros by a bandpass filter 62. with passing frequency components around a 
positive fundamental frequency of the input clock signal If necessary, the 

10 level of the extracted frequency components is doubled to compensate the 
energy of the cut off negative frequency components. An output of the 
bandpass filter 62 is inverse-transformed into the time domain using, for 
example, Inverse Fast Fourier Transform (IFFT) by the frequency domain to 
the time domain transforming part 63 to obtain an analytic signal z c (t). 

15 Furthermore, another example of the analytic signal transforming 

part 23 will be explained with reference to Fig. 14. 

A digitized input signal is stored in a buffer memory 71 . A portion of 
the signal stored in the buffer memory 71 is taken out therefrom by a signal 
talcing out part 72 in the sequ e ntial ord e r such that th e portion of th e signal 

20 b e ing tak e n out is partially ov e rlapp e d with th e portion of th e signal taken 
out just b e fore. — Th e taken out s ignal is multipli e d by a window function by 
a window function multiplying part 73, and an output signal of the window 
function multiplying part 73 is transformed into a both sided spectrum 
signal in frequ e ncy domain using Fast Fouri e r Transform by a time domain 

25 to frequency domain transforming part 1\ . — Negative frequency 

component s of this spectrum signal are made zero to obtain a single sided 
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sp o ctrum signal, and th e n r e garding this s ingl e sided spectrum signal, 
compon e nts around a fundamental fr e qu e ncy of th e input signal are retain e d 
and the oth e r fr e quency components are mad e zero by a bandwidth limiting 
part 75. — This band limited signal is transformed into a s ignal in time 
5 domain by a fr e quency domain to time domain transforming part 76 using 
Inverse FFT signal sectioning part 72 in the sequential order such that the 
current signal section and previous signal section are overlapped with each 
other. The signal section is multiplied by a window function multiplying 
part 73. and an output signal of the window function multiplying part 73 is 

10 transformed into a both-sided spectrum signal in frequency domain using 
Fast Fourier Transform by the time domain to the frequency domain 
transforming part 74. Negative frequency components of this spectrum 
signal are made zeros to obtain a single-sided spectrum signal. 
Furthermore, regarding this single-sided spectrum signal, components 

15 around a fundamental frequency of the input signal are retained and the 

other frequency components are made zeros by a bandwidth limiting part 75. 
This band-limited signal is transformed into the time domain by the 
frequency domain to the time domain transforming part 76 using Inverse 
FFT . Then this transformed signal in time domain is multiplied by an 

20 inverse window function by an amplitude correction part 77 to obtain an 
analytic signal. 

In order to decr e ase th e fr e qu e ncy of the clock signal to input a clock 
signal having a decreased frequency to the analytic signal tran s forming part 
23, as indicated by dash e d lin e s in Fig. 13, th e frequency of th e clock signal 
25 may be divided by a frequency divider 81 . — And -The input frequency range 
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of the present invention can be extended bv using a frequency divider 81. as 
illustrated bv dashed lines in Fig. 13. The fe e-frequency-divided clock 
signal may be supplied to the analytic signal transforming part 23. 
Alternatively, although not shown, using substantially jitter-free local signal, 
5 the clock signal (signal under measurement) may be converted into a signal 
with difference frequency between those signals to supply the signal to the 
analytic signal transforming part 23. 

As shown by dashed lines in Fig. 12, a low frequency component 
removing part 82 can be inserted into the output side of the linear phase 

10 removing part 27 in series so that low frequency components of the 
instantaneous phase noise A<|>(t) are removed therefrom, and the 
instantaneous phase noise A<|>(t) from which low fr e qu e ncy compon e nts are 
remov e d is supplied to the zero-crossing sampler 43. In this manner, it is 
desirable to remove low frequency components whose frequencies are 

15 sufficiently low compared with the frequency fo of the input signal, for 
example, to remove frequency components of several kHz in the case of 10 
MHz clock so that peak-to-peak jitter is not overestimated. 

Next, an embodiment of the method according to the present invention 
will be explained. Fig. 15 shows a flow-chart of the embodiment. This is, 

20 an example of the measuring method using the apparatus shown in Fig. 12. 
First in step 201, the input signal (signal under measurement) is transformed 
into a band-limited analytic signal by the analytic signal transforming part 
23. Next in step 202, an instantaneous phase of the input signal is 
estimated using the analytic signal by the instantaneous phase estimating 

25 part 26, and then in step 203, a linear instantaneous phase corresponding to 
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an ideal clock signal is estimated from this instantaneous phase by the linear 
phase estimating part 29 (Fig. 6), and the fundamental period T 0 of the 
signal under measurement is obtained by the fundamental period estimating 
part 48 from the inclination of the linear instantaneous phase. In step 204, 
5 a linear phase component is removed from the instantaneous phase by the 
linear phase removing part 27 to estimate an instantaneous phase noise A<j)(t) 
of the input signal. 

In step 205, data of the instantaneous phase noise A(|>(t) close to 
zero-crossing timings of a real part of the analytic signal are sampled by the 
10 zero-crossing sampler 43 to estimate a timing jitter sequence of the input 
signal. In step 206, a difference between approximated zero-crossing 
points from the zero-crossing point estimating part 45 is calculated by the 
zero-crossing interval calculating part 47 to estimate a zero-crossing time 
interval sequence. 

15 In step 207, a difference sequence of the timing jitter sequence is 

calculated by the differentiator 46 to estimate a period jitter sequence of the 
signal under measurement. In step 208, the period jitter sequence is 
multiplied by a ratio of the fundamental period T 0 and the zero-crossing 
time interval T^k+i by the corrector part 49 to correct the difference-based 

20 period jitter sequence. In step 209, in the state that the switch 53 is 

connected to the output side of the corrector part 49, a period jitter of the 
signal under measurement is obtained by the jitter detecting part 54 from the 
corrected period jitter sequence. 

In step 209, the peak-to-peak detecting part 32 obtains a 

25 peak-to-peak value J pp of period jitter using the equation (8), the RMS 
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detecting part 33 obtains an RMS value Jrm S of period jitter using the 
equation (7), and the 

histogram estimating part 18 obtains a histogram from the period jitter 
sequence. 

5 In step 210, in the state that the switch 53 is connected to the 

cycle-to-cycle period jitter estimating part 52, a difference sequence of the 
corrected period jitter sequence is calculated by the cycle-to-cycle period 
jitter estimating part 52 to obtain a cycle-to-cycle period jitter sequence of 
the signal under measurement. In step 211, a cycle-to-cycle period jitter of 

10 the signal under measurement is obtained by the jitter detecting part 54 from 
the cycle-to-cycle period jitter sequence. In this case, the peak-to-peak 
detecting part 32 obtains a peak-to-peak value J CC) pp of cycle-to-cycle period 
jitter using the equation (11), the RMS detecting part 33 obtains an RMS 
value Jcc,rms of cycle-to-cycle period jitter using the equation (10), and the 

15 histogram estimating part 18 obtains a histogram of the cycle-to-cycle 
period jitter. 

The estimation of the fundamental period T 0 in the step 203 and the 
calculation of the zero-crossing time interval Tk,k+i in step 206 may be 
performed before the correction in the step 208, and therefore the sequence 

20 of those processes is not limited to that of the example described above. In 
addition, the estimation of the fundamental period T 0 may be obtained 
directly from the signal under measurement. In the case where only period 
jitter is measured, the steps 210 and 211 may be omitted. In the case where 
only cycle-to-cycle period jitter is measured, the step 209 may be omitted. 

25 In the steps 209 and 211, any one or two of the peak-to-peak value, the RMS 
value, and the histogram may only be obtained. 
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In the above description, a clock signal of a microprocessor is mainly 
discussed as a signal under measurement (input signal). However, the 
present invention can be applied to a measurement of period jitter and/or 
cycle-to-cycle period jitter of a clock signal used in another equipment, or of 
5 a periodic signal such as sine wave or the like generated by another 

equipment. In addition, the input signal may be processed in the analog 
signal form rather than immediately converting it to a digital signal using 
the AD converter, and thereafter the signal may be converted, in an 
appropriate processing stage, to a digital signal. The apparatus shown in 

10 Fig. 12 may also be functioned by executing a program in a computer. 

In the above description, an instantaneous phase noise A(|>(t) is 
sampled at approximated zero-crossing points to obtain a timing jitter 
sequence A(|)[n]. However, since the liner phase remover 27 has a 
configuration shown in Fig. 6, the sampling process at the approximated 

15 zero-crossing points may be inserted in series, for example as indicated by 
dashed lines in Fig. 16, between the instantaneous phase estimating part 26 
and the continuous phase converting part 28. Alternatively, the sampling 
process at the approximated zero-crossing points may be inserted in series 
between the continuous phase converting part 28 and the linear phase 

20 estimator 29/subtractor 31. By one of those configurations, the timing 
jitter sequence A<|)[n] can also be obtained from the subtracter 31. 

In addition, the estimation of the instantaneous phase noise A(()(t) 
from an instantaneous phase is performed by the configuration of the linear 
phase remover 27 shown in Fig. 6. Therefore, the processing procedure is, 

25 as shown in Fig. 17, that after an instantaneous phase is obtained in the step 
202 in Fig. 15, the instantaneous phase is transformed to a continuous phase 
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by the continuous phase converting part 28 in step 203a, an instantaneous 
linear phase of the continuous phase is estimated by the linear phase 
estimator 29 in step 203b, and thereafter in step 204, the instantaneous 
linear phase is removed from the continuous instantaneous phase by the 
5 subtractor 31 to obtain an instantaneous noise phase A<j)(t). 

Therefore, similarly to the procedure shown in Fig. 16, as shown in 
Fig. 17, after the step 202, approximated zero-crossing samplings may be 
applied to the instantaneous phase in step 301 to obtain a sample sequence 
of instantaneous phase, and then the process moves to step 203a. In the 

10 step 203a, the sample sequence may be transformed to a continuous 
instantaneous phase. 

Alternatively, the continuous instantaneous phase obtained in the 
step 203a may be sampled at approximated zero-crossing points to obtain a 
sample sequence of continuous instantaneous phase, and then the process 

15 moves to st e p 303b step 203b . In the step 203b, an instantaneous linear 
phase may be estimated from the sample sequence of continuous 
instantaneous phase. In either case, in step 204, a timing jitter sequence 
A(|)[n] that is produced by sampling the instantaneous phase noise at 
approximated zero-crossing points. 

20 As mentioned above, according to the present invention, estimation 

errors caused by the sampling at approximated zero-crossing points can be 
decreased, a measurement result compatible with the conventional time 
interval analyzer method can be obtained, and in addition, the measurement 
can be performed in short time compared with the conventional time interval 

25 analyzer method. 
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WHAT IS CLAIMED IS: 

1 . An apparatus for measuring a jitter of a signal under measurement 
comprising: 

a timing jitter estimator to which the signal under measurement is 
inputted for obtaining its instantaneou s pha se noi s e, and to which the 
in s tantaneous phase noi s e i s inputt e d for outputting a timing jitter sequence 
that is produced by sampling the in s tantaneous phase noise at timings 
(hereinafter referred to as approximated zero crossing points) clos e to 
corresponding zero - crossing timing s of the signal under measurement timing 
jitter sequence: : 

a first differentiator to which the timing jitter sequence is inputted 
for calculating its difference sequence to output a period jitter sequence; 

a corrector part to which the period jitter sequence is inputted for 
multiplying the period jitter sequence by a ratio T 0 /Tk,k+i of a fundamental 
period T 0 of the signal under measurement and the approximated 
zero-crossing point interval Tk,k+i to output a corrected period jitter 
sequence; and 

a jitter detecting part to which the corrected period jitter sequence is 
inputted for obtaining a jitter of the signal under measurement. 

2. The apparatus for measuring a jitter of the signal under measurement 
according to claim 1 further including a second differentiator to which the 
corrected period jitter sequence is inputted for calculating its difference 
sequence to output the difference sequence to said jitter detecting part as a 
cycle-to-cycle period jitter, said second differentiator being inserted between 
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said correcting part and said jitter detecting part. 

3. The apparatus for measuring a jitter of the signal under measurement 
according to claim 1 or 2 further including a zero-crossing point detecting 
part to which a real part of the analytic signal is inputted for obtaining a 
point close to its zero-crossing timing to outpu t the sampling pul s es the 
sampling timing sequence , wherein said timing jitter estimator comprises; 

an analytic signal transforming part to which the signal under 
measurement is inputted for transforming the signal under measurement to a 
complex analytic signal; 

an instantaneous phase estimating part to which the analytic signal is 
inputted for obtaining an instantaneous phase of the analytic signal; 

a continuous phase converter for converting the instantaneous phase 
to a continuous instantaneous phase; 

a linear phase estimating part to which the continuous instantaneous 
phase is inputted for obtaining a linear phase of the continuous 
instantaneous phase; 

a subtracting part to which the linear phase and the instantaneous 
phase are inputted for removing the linear phase from the instantaneous 
phase to obtain the instantaneous phase noise; and 

a zero-crossing sampler for sampling the inputted signal using the- 
sampling puls e s t he sampling timing sequence t o output the sampled signal, 
said zero-crossing sampler being inserted in series to any one of connection 
points between said instantaneous phase estimating part and said continuous 
phase converting part, between said continuous phase converting part and 
said linear phase estimating part/subtracting part, and between said 
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subtracting part and said first differentiator. 

4. The apparatus for measuring a jitter of the signal under measurement 
according to claim 3 further including a zero-crossing interval calculating 
part to which an output timing sequence of th e sampling pul se s t he sampling 
timing sequence is inputted from said zero-crossing point detecting part for 
calculating its difference sequence to obtain the approximated zero-crossing 
point intervals T^k+i in the sequential order, and for outputting the 
approximated zero-crossing point intervals T^k+i to said corrector part. 

5. The apparatus for measuring a jitter of the signal under measurement 
according to claim 4 further including a fundamental period estimating part 
to which the linear phase is inputted from said linear phase estimating part 
for obtaining the fundamental period To from its inclination, and for 
outputting the fundamental period T 0 to said correcting part. 

6. The apparatus for measuring a jitter of the signal under 
measurement according to claim 4 further including a fundamental period 
estimating part to which the signal under measurement is inputted for 
obtaining its fundamental period To, and for outputting the fundamental 
period T 0 to said correcting part. 

7. The apparatus for measuring a jitter of the signal under measurement 
according to claim 3 further including a waveform clipper to which the 
signal under measurement is inputted for removing its amplitude modulation 
components in the state that its phase modulation components are kept 
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therein to supply the signal under measurement from which the amplitude 
modulation components have been removed to said instantaneous phase 
noise detecting part. 

8. A method of measuring a jitter of a signal under measurement 
comprising: 

a step of obtaining an instantaneous phas e noise a timing jitter 
sequence of the signal under measurement; 

a timing jitter e s timating step for s ampling the instantaneous pha s e 
noise at a timing (hereinafter referred to as approximated zero - cro s sing 
point) close to a zero cro s sing timing of the signal under mea s urem e nt to 
generate a timing jitter s e quence. ; 

a step of calculating a difference sequence of the timing jitter 
sequence to generate a period jitter sequence: 

a step of multiplying the period jitter sequence by a ratio T 0 /Tk,k+i of 
a fundamental period T 0 and the approximated zero-crossing point interval 
Tk,k+i to obtain a corrected period jitter sequence; and 

a step of obtaining a period jitter of the signal under measurement 
from the corrected period jitter sequence. 

9. A method of measuring a jitter of a signal under measurement 
comprising: 

a step of obtaining an in s tantaneou s phase noise a timing jitter 
sequence of the signal under measurement; 

a timing jitter e stimating stop for sampling th e instantaneou s phase 
nois e at a timing (h e r e inaft e r referred to as approximated z e ro crossing 
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point) clos e to a zero cro ss ing timing of th e signal under m e asur e ment to 
g e n e rate a timing jitter sequence; 

a step of calculating a difference sequence of the timing jitter 
sequence to generate a period jitter sequence: 

a step of multiplying the period jitter sequence by a ratio To/Tk,k+i of 
a fundamental period To and the approximated zero-crossing point interval 
Tk,k+i to obtain a corrected period jitter sequence; 

a step of calculating a difference sequence of the corrected period 
jitter sequence to generate a cycle-to-cycle period jitter sequence: and 

a step of obtaining a cycle-to-cycle period jitter of the signal under 
measurement from the cycle-to-cycle period jitter sequence. 

10. The method of measuring a jitter of the signal under measurement 
according to claim 8 or 9 further including a step of obtaining a point close 
to a zero-crossing timing of a real part of the analytic signal to obtain the 
approximated zero-crossing point, wherein said timing jitter estimating step 
comprises: 

a step of transforming the signal under measurement to a complex 
analytic signal; 

a step of obtaining an instantaneous phase of the signal under 
measurement from the analytic signal; 

a step of transforming the instantaneous phase to a continuous 
instantaneous phase; 

a step of obtaining a linear phase from the continuous instantaneous 

phase; 

a step of removing the linear phase from the continuous 
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instantaneous phase to obtain the instantaneous phase noise; and 

a step of sampling any one of the instantaneous phase, the continuous 
instantaneous phase and the instantaneous phase noise at the approximated 
zero-crossing timing. 

1 1 . The method of measuring a jitter of the signal under measurement 
according to claim 10 further including a step of calculating a difference 
sequence of a timing sequence that represents the approximated 
zero-crossing points to obtain the approximated zero-crossing point intervals 
T k k+i in the sequential order. 

12. The method of measuring a jitter of the signal under measurement 
according to claim 1 1 further including a step of obtaining the fundamental 
period T 0 from an inclination of the linear phase. 

13. The method of measuring a jitter of the signal under measurement 
according to claim 1 1 further including a step of obtaining the fundamental 
period To from the signal under measurement. 

14. The method of measuring a jitter of the signal under measurement 
according to claim 10 further including a step of removing amplitude 
modulation components of the signal under measurement in the state that its 
phase modulation components are kept therein to move to said timing jitter 
estimating step. 
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ABSTRACT 

A signal under measurement x(t) is transformed into a complex 
analytic signal z c (t), and an instantaneous phase of the x c (t) is estimated 
using the z c (t). A linear phase is removed from the instantaneous phase to 
obtain a phase noise waveform A(f>(t) of the x(t), and the A(j)(t) is sampled at 
a timing close to a zero-crossing timing of the x(t) to obtain a timing jitter 
sequence. Then a difference sequence of the timing jitter sequence is 
calculated to obtain a period jitter sequence. The period jitter sequence is 
multiplied by a ratio T 0 /Tk,k+i of the fundamental period T 0 of the x(t) and 
the sampling time interval Tk,k+i to make a correction of the period jitter 
sequence. A period jitter value of the x(t) is obtained from the corrected 
period jitter sequence. 



